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Proximity to a Mott insulating phase is likely to be an important physical ingredient of a the- 
ory that aims to describe high-temperature superconductivity in the cuprates. Quantum cluster 
methods are well suited to describe the Mott phase. Hence, as a step towards a quantitative theory 
of the competition between antiferromagnetism and d-wave superconductivity in the cuprates, we 
use Cellular Dynamical Mean Field Theory to compute zero temperature properties of the two- 
dimensional square lattice Hubbard model. The d-wave order parameter is found to scale like the 
superexchange coupling J for on-site interaction U comparable to or larger than the bandwidth. 
The order parameter also assumes a dome shape as a function of doping while, by contrast, the 
gap in the single-particle density of states decreases monotonically with increasing doping. In the 
presence of a finite second-neighbor hopping t' , the zero temperature phase diagram displays the 
electron-hole asymmetric competition between antiferromagnetism and superconductivity that is 
observed experimentally in the cuprates. Adding realistic third-neighbor hopping t" improves the 
overall agreement with the experimental phase diagram. Since band parameters can vary depending 
on the specific cuprate considered, the sensitivity of the theoretical phase diagram to band param- 
eters challenges the commonly held assumption that the doping vs Tc/TJJ^'^^ phase diagram of the 
cuprates is universal. The calculated angle-resolved photoemission spectrum displays the observed 
electron-hole asymmetry. The tendency to homogeneous coexistence of the superconducting and an- 
tiferromagnetic order parameters is stronger than observed in most experiments but consistent with 
many theoretical results and with experiments in some layered high-temperature superconductors. 
Clearly, our calculations reproduce important features of d-wave superconductivity in the cuprates 
that would otherwise be considered anomalous from the point of view of the standard Bardeen- 
Cooper-Schrieff'er approach. At strong coupling, d-wave superconductivity and antiferromagnetism 
appear naturally as two equally important competing instabilities of the normal phase of the same 
underlying Hamiltonian. 

PACS numbers: 71.10.-w, 71.35.-y, 72.80.Sk, 78.30.Am 



I. INTRODUCTION 



Superconductivity in the cuprates and in the layered 
organics of the BEDT family is highly anomalous, i.e., it 
displays a number of properties that cannot be explained 
by the Bardeen-Cooper-Schrieffer (BCS) theory modified 
for d-wave symmetry. For example in the cuprates, su- 
perconductivity emerges upon doping an antiferromag- 
netic Mott insulator. Moreover, in the so-called under- 
doped region near the insulator, experiments show that 
the gap in the single-particle density of states decreases 
upon doping while Tc/T™°^, or the order parameter, in- 
creases, in sharp contrast with expectations from stan- 
dard BCS theorjii. In the organics, antiferromagnetism 
and superconductivity are separated by a first order tran- 
sition and a Mott transition separates the correspond- 
ing states with no order. For both the organics and the 
cuprates, there is much evidence from approximate solu- 
tions that the essential low-energy physics is described by 
the one-band Hubbard model for the appropriate lattice, 
band-structure, interaction and dopings^^i^. 



Understanding anomalous superconductivity theoreti- 
cally in a quantitative manner is still a challenge. An 
important step towards this goal is to obtain accurate 
solutions of the Hubbard model. Despite the apparent 
simplicity of the model, it is extremely difficult to solve 
in the relevant regime where neither potential [U) nor 
kinetic energy (8t) dominate. In recent years, a number 
of numerical methods have shed light on this problem. In 
this paper, we will describe the results obtained from Cel- 
lular Dynamical Mean-Field Theorj*^ (CDMFT) for d- 
wave superconductivity and its competition with antifer- 
romagnetism in the cuprates. The corresponding study 
in the organics has been published^ii. The results will 
be compared with other quantum cluster methods^i^i^, 
mean-field theories, slave boson, and with variational ap- 
proaches. 

Our choice of method is motivated by the following 
considerations. Since anomalous superconductivity ap- 
pears near antiferromagnetic Mott insulating phases, it 
is important to use approaches that treat these phases 
correctly. CDMFT is a generalization of Dynamical 
Mean-Field Theory (DMFT>W'i2_ ^j^g jitter method de- 



2 



scribes the Mott insulator-metal transition of the Hub- 
bard model exactly in the limit of infinite dimensions. 
DMFT is, by construction, a local theory that maps the 
full interacting many-body lattice problem onto a single- 
site embedded in a self-consistent bath. Unfortunately, 
the local nature of the spatial correlations inherent to 
single-site DMFT precludes a description of the super- 
conducting phase with c?-wave symmetry observed ex- 
perimentally in the cuprates. This limitation has been 
overcome by several recently developed cluster extensions 
of DMFT. In addition to CDMFT^ these include the 
Dynamical Cluster Approximation (DCA)^^ and Varia- 
tional Cluster Perturbation Theory (VCPT)^'!^^ also 
known as Variational Cluster Approximation (VGA). 
These methods incorporate short-range correlations as 
well as some properties of the infinite lattice in a system- 
atic and causal manner (for reviews see^i^iiiSi) . CDMFT 
maps a lattice problem onto a finite-size cluster (with 
different boundary conditions than DCA) and is able to 
describe the short-range correlations accurately within 
the cluster. A self-consistently determined bath of un- 
correlated electrons approximates the effect of the rest of 
the infinite lattice on the cluster. Coupling between the 
embedded cluster and the bath provides a self-consistent 
theory that naturally allows for phase transitions and 
phases with long range order. In a way, CDMFT makes 
a compromise between short-range and long-range corre- 
lations, the latter being associated with order parameters 
that appear only in the bath. 

Cluster models can be formulated in a general func- 
tional frameworkiSii^ and they all converge to the same 
limit for very large cluster sizes. Nevertheless, different 
cluster methods have convergence rates that depend on 
the physical observable. In particular, at finite tempera- 
tures, CDMFT has been shown to converge exponentially 
quickly for local quantities such as the density of statesi^. 
The quality of the approximation made by CDMFT has 
been tested extensively by comparisons with exact re- 
sultsiiii^. The method was found to be accurate in de- 
scribing the complexities of the Mott insulator-to-metal 
transition, indicating that both high as well as low energy 
phenomena are well captured. Furthermore, CDMFT 
has also been used successfully to elucidate the differ- 
ences in the evolution of a Mott-Hubbard insulator into 
an anomalous correlated metal for the electron- and hole- 
doped casesi2i2£ as well as to shed light on the normal 
state pseudogap phenomenon3i. 

In the following section, we describe our implemen- 
tation of CDMFT in more details. We then present 
the zero-temperature results and compare with those of 
other numerical approaches. It is shown that, in the in- 
termediate coupling regime U — 8t, the d-wave super- 
conducting phase at zero temperature is stabilized. In 
addition, at strong coupling the magnitude of the order 
parameter scales with J = At^ /U. We also discuss the ef- 
fect of longer-range hopping and of the proximity to the 
Mott transition. We present the phase diagram for both 
hole and electron-doped systems, exhibiting the compe- 



tition between d-wave superconductivity and antiferro- 
magnetism. The regions where d-wave superconductivity 
and antiferromagnetism appear are in semi-quantitative 
agreement with experiment, except for a region of coex- 
istence that is not always observed experimentally. It 
is striking to observe that as we make the model more 
realistic by including all the main band structure param- 
eters t, t' and , the coexistence region decreases and the 
overall phase diagram becomes closer to experiments. Fi- 
nally, before we summarize the results, we compute the 
single-particle density of states and show that the or- 
der parameter increases with doping in contrast with the 
single-particle gap that decreases, in qualitative agree- 
ment with experiment. In addition, we find that the 
results of Angle Resolved Photoemission Spectroscopy 
(ARPES)'^^ have a natural explanation in our approach. 
Discussion of our results in the context of a sample of the 
existing literature appears at the end of each subsection. 



II. MODEL AND METHODS 



We define the Hubbard model Hamiltonian by 



H 



Uj2n 



(1) 



where and U correspond to the hopping and the on- 
site Coulomb repulsion respectively. In order to investi- 
gate the d-wave superconducting phase within CDMFT 
we write an effective action containing a Weiss dynamical 
field Qo with both normal (particle-hole) and anomalous 
(particle-particle) components that describes the degrees 
of freedom outside the cluster (the bath) as a time de- 
pendent hopping within the cluster 
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For the case of a 2 x 2 plaquette, which we shall con- 
sider throughout this work, the Nambu spinor is defind 
by ^fj; = (fij|, . . . , ^41, . . . , d4x), and the greek letters 
/i, V label the degrees of freedom within the cluster. Phys- 
ically, this action corresponds to a cluster embedded in a 
self-consistently determined medium with c?-wave pairing 
correlations. 

Given the effective action with a starting guess for the 
Weiss field we compute the cluster propagator Gc by 
solving the cluster impurity Hamiltonian that will be de- 
scribed shorty. Then we extract the cluster self energy 
from tc = Qq^ - G-\ Here, 

is an 8 X 8 matrix, G^^^^ = -{Tdfj_cr{T)dl^{0)) and = 
— {Tdf^l(T)d^l (0)) are the imaginary-time ordered normal 
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and anomalous Green functions respectively. Using the 
self-consistency condition, 



with 



Nr. 



{2n) 



c?k G(k, iujn) 



+ Se(«C^„) (4) 



(5) 



we recompute the Weiss field Go and iterate till conver- 
gence. Here i(k) is the Fourier transform of the super- 
lattice hopping matrix with appropriate sign flip between 
propagators for up and down spin and the integral over 
k is performed over the reduced Brillouin zone of the 
superlattice. 

To solve the cluster impurity problem represented by 
the effective action above, we express it in the form of a 
Hamiltonian -ffimp with a discrete number of bath orbitals 
coupled to the cluster and use the exact diagonalization 
technique (Lanczos method)^ 



p.va ma 

+ ^mM<Tam<T(cA<<7 + h.C.) + C/ ^ "mT'^mX 

a 

-|-a2|a"| — ag^ajj + a'i^a'^i — + h.c). 



Here H^v = l,...,Nc label the sites in the cluster and 
Efiva represents the hopping and the chemical potential 
within the cluster. The energy levels in the bath are 
grouped into replicas of the cluster {Nc = 4) (two replicas 
in the present case) with the labels m = 1, - ■ ■ , Nc and 
a = 1,2 such that we have 16 bath energy levels e^^. 
coupled to the cluster via the bath-cluster hybridization 
matrix V^^^. Using lattice symmetries we take V^^^ = 
V^Sm/j, and e^^ = e". The quantity A" represents the 
amplitude of superconducting correlations in the bath. 
No static mean-field order parameter acts directly on the 
cluster sites^. 

The parameters e", V" and A" are determined by im- 
posing the self-consistency condition in Eq.|4]using a con- 
jugate gradient minimization algorithm with a distance 
function 



Y (p'o ^ i'i-^n) ~ Go ^ (iuJn) 



fiiy 



(6) 



that emphasizes the lowest frequencies of the Weiss field 
by imposing a sharp cutoff at ojc = 1.5. (Energies are 
given in units of hopping t, and we take h = 1 and = 
1.) The distance function in Eq.® is computed on the 
imaginary frequency axis (effective inverse temperature, 
P = 50) since the Weiss field Go{'ii^7i) is a smooth function 
on that axis. 



With the bond superconducting order parameter de- 
fined as 



(7) 



we consider c?-wave singlet pairing {tp = -012 = —'023 — 
■034 — — "041). The average is taken in the ground state 
of the cluster. 

To study the competition with antiferromagnetism, 
while preserving bipartite symmetry, hybridization 
V^^^ = V^Smfi and bath site energies e^^ = e" become 
spin dependent. This doubles the number of indepen- 
dent hybridization and bath parameters. The staggered 
magnetic order parameter is given by the cluster wave 
function average 



(8) 



with M = Ml ^ -Ah = Ah = -M4. 

The finite size of the bath in the exact-diagonalization 
technique is an additional approximation to the CDMFT 
scheme. The accuracy of this approximation can be ver- 
ified by comparing the CDMFT solution for the one- 
band Hubbard model with the solution from the Bethe 
ansatzi^i^. We have also used this comparison in one 
dimension as a guideline to fix the choice of parameters 
in the distance function {ojc = 1-5 and (3 = 50). These 
results in one dimension also compare well with those 
obtained using the Hirsch-Fye Quantum Monte Carlo 
algorithm as an impurity solver where the bath is not 
truncated2£. Further, using finite-size scaling for these 
low (but finite) temperature calculations^^, it was shown 
that, at intermediate to strong coupling, a 2 x 2 cluster 
in a bath accounts for more than 95% of the correla- 
tion effect of the infinite size cluster in the single-particle 
spectrum. 

We can also perform an internal consistency check on 
the effect of the finite bath on the accuracy of the calcu- 
lation. With an infinite bath, convergence insures that 
the density inside the cluster is identical to the den- 
sity computed from the lattice Green function. In prac- 
tice, we find that there can be a difference of ±0.02 be- 
tween the density estimated from the lattice and that 
estimated from the cluster. We will display results as 
a function of cluster density since benchmarks with the 
one-dimensional Hubbard model show that, with a finite 
bath and the procedure described above, one can repro- 
duce quite accurately Bethe ansatz results for n (/i) when 
the cluster density is used. Nevertheless, we should adopt 
a conservative attitude and keep in mind the error esti- 
mate mentioned above. 



III. ANOMALOUS SUPERCONDUCTIVITY 
AND ITS COMPETITION WITH MAGNETISM 

In this section we first discuss d-wave superconductiv- 
ity by itself. Competition with magnetism is taken into 
account in subsequent subsections. We will conclude with 
single-particle spectral properties. 
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FIG. 1: (Color online) a) d-wave order parameter i/) as a func- 
tion of filling n and onsite Coulomb repulsion U, t' = 0. b) 
Same results but with vertical axis divided by J = 4t^ /U. 



A. Scaling as J at strong coupling 

The effect of interaction strengtli on d-wave supercon- 
ductivity is illustrated in Fig. [T] a, where we show the 
d-wavc order parameter ip as a function of the density 
n for t' — and different values of onsite Coulomb re- 
pulsion U. No long-range antiferromagnetism is allowed 
in this calculation. As U is increased to 8t, the or- 
der parameter acquires a very broad maximum around 
S = |1 — n| = 0.15 (optimal doping). The c?-wave Weiss 
field (not shown), on the other hand, is monotonic as a 
function of doping with a maximum value close to half- 
filling. As the system approaches half-filling from opti- 
mal doping there is a progressive localization due to the 
increasing proximity of the Mott insulator and hence, de- 
spite the stronger c?-wave Weiss field in this region, there 
is a suppression of the d-wave order parameter. This 
suppression is not seen for U — At that is too small for 
a Mott transition. However, even in that case, if we 
were to allow for long-range antiferromagnetic correla- 
tions, the insulating gap at half-filling that emerges from 
antiferromagnetic (Slater) correlations should also sup- 
press d-wave superconductivity completel y^^'^^ . Clearly, 
antiferromagnetism near half-filling can destroy super- 
conductivity but the paramagnetic Mott transition by 
itself, without antiferromagnetism, suffices at strong to 
intermediate coupling. The maximum of the order pa- 
rameter increases as U increases but then it decreases as 
U changes to \2t, l&t and 24i. This clearly signals that d- 
wave superconductivity will be strongest at intermediate 
coupling and that at strong coupling the d-wave order pa- 
rameter scales as the superexchange coupling J = At'^/U. 
In Fig. [1] b, the value of the order parameter is divided 
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FIG. 2: d-wave order parameter for [/ = 8t as a function of 
filling n for various values of t' /t (frustration). 



by J on the vertical axis to better exhibit this scaling 
for U > I2t. The value of optimal doping is nearly in- 
dependent of U in the intermediate to strong coupling 
regime. 

On the overdoped side, d-wave superconductivity dis- 
appears at a doping that is comparable with the Dy- 
namical Cluster Approximation (DCA)^^ and YCAM. 
Variational calculations^!^ find d-wave superconductiv- 
ity up to about 35% doping. The strong-coupling scaling 
with the superexchange coupling J = At'^/U as well as 
the dome like shape of the superconducting order param- 
eter, which are a hallmark of correlated superconductiv- 
ity, were also found in VCA5S,, renormalized mean-field 
theorji^^, gauge theories^ and in early RVB slave bo- 
son studies^i^^. In multiorbital models for fuUerenes, a 
dome behavior is also found as a function of the corre- 
lation strengt h^ ''i^^ . We note that at large doping for 
U = 8t— 12t, we find in Fig. [T] that the fall of the d-wave 
order parameter occurs with constant negative curvature 
as a function of doping, as seems to occur experimentally. 
This contrasts with weak coupling studie s^'''"?- where this 
does not occur. 

Recent experiments^ in four families (dif- 
ferent x) of high-temperature superconductors 
(Caa;Lai_^)(Bai.75_2;Lao.25+2:)Cu30y suggest that 
scaling with J has been observed experimentally in the 
cuprates. 



B. Effect of second- neighbor hopping t' /t without 
competition with antiferromagnetism 

The effect of next- nearest- neighbor hopping t' /t is il- 
lustrated in Fig. [5] where we plot the d-wave order param- 
eter ^ as a function of doping for various values of the 
next-nearest neighbor (diagonal) hopping t' . For clarity, 
we represent both the hole and electron-doped cases on 
the same plot by performing a particle-hole transforma- 
tion and considering t' — +0.3t in the electron-doped 
case. The maximum of ip grows with increasing t'. 

This trend with t' was found in DCA^^ , in recent vari- 
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ational studies^ and in two-leg ladder studies of the 
t — t' — J model^. However, the increase of -0 with 
t' does not seem consistent with the empirical correla- 
tion found in Ref. between a larger optimal Tc and a 
smaller value of the t' obtained from band structure cal- 
culations. As found in an earlier VGA study^ the maxi- 
mum ip is larger for electron doping than for hole doping 
when calculations are done on 2 x 2 clusters. However, 
the situation reverses for larger clusters. This suggests 
that finite-size effects may be influencing our conclusions 
on the effect of t' . Superconductivity is known to be 
strongest on the hole- as opposed to the electron-doped 
compounds (Using particle-hole transformation, the lat- 
ter correspond to positive t' on Fig (2]). Nevertheless, we 
will see in the following section that including the com- 
petition with antiferromagnetism restores the experimen- 
tally observed electron-hole asymmetry. 



C. Phase diagram and competition with magnetism 

In the previous sections, we have constrained the bath 
of the associated impurity problem to allow only normal 
state and d-wave superconducting solutions. In this sec- 
tion, we start from initial bath parameters that allow for 
both d-wave and antiferromagnetic symmetry breaking. 
Fig. [3] displays the results for the zero-temperature phase 
diagram for t' = — 0.3t and three values of the interac- 
tion strength U — 6t, U — 8t and U — 12t. The plots are 
as a function of electron density, so that the system is 
hole-doped for n < 1 and electron-doped for n > 1. The 
order parameter for d-wave superconductivity is multi- 
plied by a factor of 10 so that it is on the same scale as 
the antiferromagnetic order parameter. We have checked 
that in the coexistence region, even if we start the iter- 
ations with a fully converged antiferromagnetic solution, 
the final solution is the same as the one exhibited in Fig. 

m 

From Fig. [31 we observe that antiferromagnetism oc- 
curs over a narrower range of dopings as U increases 
since J then decreases (the trend would be opposite at 
weak coupling). There is homogeneous coexistence of 
antiferromagnetism and d-wave superconductivity near 
half-filling in all cases. That phase can be called a su- 
perconducting antiferromagnet or an antiferromagnetic 
superconductor. D-waye superconductivity exists by it- 
self at large electron or hole doping. The transition from 
a homogeneous coexistence phase to pure d-wave super- 
conductivity is second order, except for large values of 
U = 12t on the electron-doped side. Qualitatively, it 
seems that U ~ 8t gives a better agreement with the 
experimental phase diagram of the cuprates since, in 
that case, superconductivity appears alone over a broader 
range of dopings on the hole than on the electron-doped 
side. Also, the value of the maximum d-wave order pa- 
rameter is larger on the hole- than on the electron-doped 
side, showing that competition with antiferromagnetism 
can reverse the trend observed as a function of t' in Figl^l 




FIG. 3: (Color online) Superconducting (red circles) and anti- 
ferromagnetic (blue squares) order parameters for t' = — 0.3t 
and, from top to bottom, the three values U = 6t, U — 8t, 
and U = 12t. The amplitude of the superconducting order 
parameter is multiplied by a factor of 10 to be on a scale 
comparable with the antiferromagnetic one. 



Choosing a value of U on the electron-doped side^i^ 
that is smaller than the value of U for the hole-doped 
side would also help in making the tendency for d-wave 
superconductivity smaller on the electron-doped side, as 
seems to be observed experimentally. The asymmetry in 
the maximum value of the d-wave order parameter for 
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hole and electron doping is also observed for U = 6t, 
but in that case the antiferromagnetism occurs over a 
doping range that is unreasonably large compared with 
experiment. For U — 8t, optimal doping occurs around 
15% in the hole-doped case. U — 8t is also consistent 
with the value necessary to explain details of the spin 
wave spectrum obtained by neutron measurements at 
half-filling^i^. All these qualitative trends agree with 
the experimental phase diagram except for the follow- 
ing: Antiferromagnetism extends over a broader range 
of dopings on the hole-doped side than observed exper- 
imentally (see however the next section) and there is a 
strong tendency for homogeneous coexistence of the two 
order parameters, even though the d-wave order parame- 
ter is suppressed by the presence of antiferromagnetic or- 
der. While the suppression of d-wave superconductivity 
by antiferromagnetism is appreciable, the reverse effect 
is almost negligible"*^. This is also observed in mean- field 
studies^. 

Homogeneous coexistence of antiferromagnetism and 
d-wave superconductivity is not generic in the cuprates. 
Nevertheless it has been observed recently in ordered 
layered compounds^ and in the electron-doped cuprate 
PrCeCuO^. Although such coexistence is not observed 
in compounds like YBCO, it appears in La2CuO4.11 at 
zero field and in Lai.9Sro.iCu04 under applied mag- 
netic field^-. The agreement of the latter field-dependent 
experiments-- with the theoretical predictions^'^ re- 
veals the proximity of single-phase d-wave superconduc- 
tivity with homogeneous coexistence of antiferromag- 
netism and d-wave superconductivity. Other LSCO 
compounds^i^ reveal the proximity of antiferromag- 
netism and d-wave superconductivity through the appli- 
cation of a magnetic field. Muon spin rotation studies as 
a function of doping, temperature and field show that 
zero-field coexistence of superconductivity with short- 
range static magnetisni^ is generic in the underdoped 
regime. Since competition between antiferromagnetism 
and d-wave superconductivity seems to be quite sensi- 
tive to disorder^, the homogeneous coexistence between 
the two phases that we find here may be reconcilable with 
experiments. 

It is instructive to compare our theoretical results with 
other approaches that include competition of d-wave su- 
perconductivity with antiferromagnetism. Early varia- 
tional calculations^ aX t' = 0, U = 10 and with the 
t — J model^ showed a strong tendency for antifer- 
romagnetism and superconductivity to coexist homoge- 
neously, even though there is a region where d-wave su- 
perconductivity exists by itself. The most recent varia- 
tional calculations^!^ show the same trend. Weak cou- 
pling calculations with functional renormalization cou- 
pled to renormalized mean- field theory at U — 2.5 and 
t' — —0.15 show also a coexistence region2§,. A re- 
cent review^ of variational approaches and renormalized 
mean-field theory results^ shows that a coexistence re- 
gion is often obtained in these approaches, as it was in 
early slave-boson theories^. 



An early version of cluster DMFTS? using t' = — 0.15t 
and U = 4.8t found that d-wave superconductivity and 
antiferromagnetism coexist over most of the doping range 
1 - n < 0.3. Recent CDMFT calculationsS at i' = 
find that an homogeneous coexistence region appears at 
weak coupling but is replaced by a first order transition 
between both phases at strong coupling. Phase separa- 
tion as a function of doping takes place in strong coupling 
but not in weak coupling. Here, at finite t', we did find 
a first-order transition at large U on the electron-doped 
side but phase coexistence persists at low doping. 

In VGA calculations that include cluster chemical po- 
tential as a variational paramete r^'^i^^ to insure thermo- 
dynamic consistency, a coexistence phase with continu- 
ous transition to d-wave superconductivity is found for 
t' = —0.3 at weak coupling, but at strong coupling the 
transition between homogeneous coexistence and pure d- 
wave superconductivity is first order^^, as we have found 
here for U = 12t on the electron-doped side. In early 
VGA studies^, U = 8, t' = -0.3, it was noticed that the 
energy scale for that first order transition is larger on the 
hole than on the electron doped side. In VGA calcula- 
tions that do not include chemical potential as a vari- 
ational parameter, zero-temperature finite-size studies 
find antiferromagnetism and d-wave superconductivity 
over doping ranges that agree semi-quantitatively with 
experiment for both hole and electron doped systems^S., 
using U = 8t and taking second nearest neighbor hopping 
parameter t' = — 0.3t and third-nearest neighbor hopping 
t" — Q.2t from band structure. As in our case, there is a 
coexistence region on the electron-doped side, but on the 
hole-doped side there is still some non-monotonic size de- 
pendence that suggests that inhomogeneous phases may 
be more stable. 

Finally, in DGA, system sizes up to 32 sites have been 
studied for U = 4t, t' = 0. At 10% doping, a super- 
conducting transition is found^* with critical tempera- 
ture Tc = 0.023. The question of ground state coex- 
istence was not studied since these are finite tempera- 
ture studies. However, it was established at intermediate 
coupling that the frequency and momentum dependence 
of the pairing channel comes from the antiferromagnetic 
fluctuations^iZOiZi. Within DGA, there is a tendency 
for phase separation in the normal state on the electron- 
doped side of the phase diagram^^. The same tendency 
is observed at weak coupling2^ using the Two-Particle 
Self-Gonsistent (TPSG) approacbZi. 



D. Additional effects of tiie band structure on the 
phase diagram 

The one-band Hubbard model is an effective model of 
the cuprates, so parameters can vary with doping^i^. 
Given our numerical uncertainties, it would be inappro- 
priate at this point to fine-tune our parameters to try 
to find a best fit of the experimental phase diagram of 
the high-temperature superconductors. Nevertheless, we 
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FIG. 4: (Color online) Superconducting (red circles) and an- 
tiferromagnetic (blue squares) order parameters for U = St 
and various values of the band parameters. On the top, pa- 
rameters t' — — 0.3t, t" — 0.2t appropriate for YBCO. On the 
bottom, t' = -0.17t, t" = 0.08t as in LSCO. The amplitude of 
the superconducting order parameter is multiplied by a factor 
of 10 to be on a comparable scale with the antiferromagnetic 
one. 



point out in this section that when more realistic band 
parameters are used for the calculation, the qualitative 
changes of the resulting phase diagram go in the right di- 
rection to improve agreement with experiment. We also 
comment on the implications of sensitivity to band pa- 
rameters. 

The most common parameters used for the one-band 
model of YBCO for examp\eB21 are f = -0.3t, t" = 
0.2t, where t" stands for third nearest neighbor hopping. 
It has been known for a decade^^ that for the half-filled 
parent compound Sr2Cu02Cl2, for example, these val- 
ues describe well ARPES experiments. Adding third- 
nearest neighbor hopping t" — 0.2t to the calculation 
of the previous section leads to the phase diagram illus- 
trated on the upper part of Fig. SI Clearly the phase dia- 
gram becomes qualitatively closer to that of the cuprates 
since the antiferromagnetic phase has a stronger electron- 
hole asymmetry and the tendency for coexistence be- 
tween antiferromagnetism and c?-wave superconductivity 



is reduced. These two results were found first in recent 
VCA calculations^i. Note also that the phase bound- 
aries for antiferromagnetism are extremely close to those 
obtained from VCA for the same parameters^. In our 
CDMFT calculations, the effect of t" enters only in the 
self-consistency condition Eq.Q through the parameter 
t(k) of the superlattice Green function Eq.®. In VCA 
the clusters are larger so that t" can also enter in the 
cluster calculation, not only in the self-consistency con- 
dition. The excellent agreement between both methods 
for the boundaries of the antiferromagnetic phases sug- 
gests that the main effect of t" is taken into account in 
our CDMFT calculations. 

On the lower part of Fig. 21 one finds the 
phase diagram corresponding to parameters typical of 
La2-:,Sr^Cu04 (LSCO)i^, namely t' = -0.17t and t" = 
O.OSt. There are quantitative changes compared with the 
top figure appropriate for YBa2Cu307_a; (YBCO). In 
particular, on the hole-doped side of the phase diagram 
the d-wave phase extends to larger dopings for the above 
LSCO parameters than those for YBCO. Although there 
is no qualitative change, this quantitative sensitivity of 
the phase diagram to band parameters suggests that the 
universality of the Tc {S) /T™"^ relatio n ''^i''^ , usually as- 
sumed for hole-doped (S) high-temperature superconduc- 
tors, may not be completely accurate, unless it signals 
that band parameters are in fact quite close, as has been 
suggested in some recent studies^. The possible lack of 
universality of the Tc (S) /T'™'*^ relation has been noticed 
before^^. In experiments, it has been pointed out a while 
ago^^ that single layer Bi2Sr2-2La2Cu206+5 (BSLCO) 
deviates from the universal Tc (S) / J"™"^ curve. The elec- 
tronic structure of that compound differ^^ from that of 
LSCO, the compound used to determine the universal 
Tc (S) /TJ""^ curved. Hence, band structure could be 
one of the factors explaining the deviations of single layer 
BSLCO. There is also a "universal" Tc{S {290 K))/T^^'''' 
determined from thermopower S{290K) at 290K instead 
of doping^-. However, thermopower has been experimen- 
tally linked to band structure^. 



E. Gap in total density of states and ARPES 
spectrum 

Having established that the d-wave order parameter '0 
decreases towards zero as we approach half-filling, we now 
show that the gap A computed from the single-particle 
density of states increases monotonically towards half- 
filling n = 1 in the hole-doped case. All the calculations 
of this section are for t' = —0.3t and U — 8t. The local 
density of states can be obtained from the cluster or from 
the lattice Green function. The results are identical in 
the limit of infinite bath size. The single-particle gap is 
extracted from the local density of states, that we obtain 
from the lattice Green function that will be used later 
to obtain the single-particle spectral weight A (k, w = 0). 
The lattice Green function is computed by periodizing 
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FIG. 5: The gap as a function of filling, for U = 8i, t' = -O.St. 
The gap is defined as half the distance between the two peaks 
on either side of lj = 0, as they appear, for example, in the 
inset. 

the superlattice Green function G^^CkjUs) in Eq.([5]) over 
the entire Brillouin zone of the original infinite lattice as 
is done in CPT— . This proceeds as follows. In reciprocal 
space, any wave vector k in the Brillouin zone may be 
written as k = k + K where both k and k are continuous 
in the infinite size limit. However k is defined only in 
the reduced Brillouin zone that corresponds to the su- 
perlattice. On the other hand, K is discrete and denotes 
reciprocal lattice vectors of the superlattice. Then the 
lattice Green function is defined by 

Giatt(k,c^) = ^e-^''-('>-'--)G''''(k,c^), (9) 

where G^'^{]i,Lu) = G^''(k,a;). Integration over k yields 
the local density of states. Different periodization proce- 
dures can be used-i^. 

The single-particle gap is defined as half the distance 
between the peaks that surround the minimum near 
a; = 0. The inset of Fig. [5] shows the local density of 
states in the superconducting state at n = 0.9. Due to 
the finite value of the broadening parameter (O.lt), some 
spectral weight fills the d-wave gap, which vanishes lin- 
early near uj — 0. We do not focus on this oj ^ Q 
energy scale. In Fig. [S] we exhibit the monotonic in- 
crease of the single-particle gap as filling increases toward 
n = I (unlike the order parameter decrease) in accor- 
dance with, for example, thermal conductivity^ and tun- 
neling experiment s^^ ' for the cuprates. This is a radi- 
cal departure from a conventional superconductor where 
the single-particle gap A and tJj are proportional to each 
other. Experimentally in the hole-doped cuprates^, the 
magnitude of the gap at optimal doping is estimated to 
be around bOmeV. Here we obtain 0.2t, which is in good 
agreement with experimental estimates if we take a rea- 
sonable value t = 250meV. CDMFT should be reHable 
for the value of this gap since it is extracted from a local 
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FIG. 6: (Color online) Hole-doped d-wave superconducting 
phase for n — 0.94, U = St, t' = —0.3. Upper figure, color 
contour plot of A (k, a; = 0^) in first quadrant of the Brillouin 
zone. Energy resolution is r; = O.lt for both figures. Lower 
figure, A(k,u!) for various values of wave vector along high- 
symmetry directions. The relative scale in the contour plot 
can be made quantitative by comparison with the lower figure. 



quantity. 

In Figs, inland [71 we present single-particle spectra in 
the d-wave superconducting state for, respectively, hole 
doping n = 0.94 and electron doping n = 1.06. The upper 
part of the figures illustrates the momentum distribution 
curves (MDC) at the Fermi energy in the first quadrant 
of the Brillouin zone, in other words a contour plot of the 
single-particle spectral weight ^ (k, w = 0) with a finite 
energy resolution rj — O.lt. The lower part contains en- 
ergy distribution curves (EDC) A (k, uj) as a function of 
frequency in units of t for various wave vectors k along 
symmetry directions. 

The MDC for hole doping, n = 0.94, in Fig. |6land elec- 
tron doping, n — 1.06, in Fig. [71 are very similar to the 
corresponding results in the normal state obtained earlier 
with the same approach (Figs. 5 and 6 of Ref. [2l|, and 
Fig. 4 of Ref. Ho) and with CPT-^^-SS^ the latter having 
the best momentum resolution. They display the striking 
asymmetry of the hole- and electron-doped cases that can 
be measured by ARPES in the c?-wave superconducting 
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FIG. 7: (Color online) Same plots as in Fig. |6] but for the 
electron-doped case n = 1.06. The maximum on the color 
scale corresponds to 104% of the maximum of the data to 
stress finer details near the maximum. All other parameters 
are identical. 

phase22 and in the pseudogap phase. More specifically, 
for the hole-doped case the Fermi energy plot (MDC) 
shows weight near (7r/2,7r/2) with a clear suppression 
everywhere else. On the other hand, the electron-doped 
case Fig. [7]shows weight in complementary areas of (tt, 0) 
and (0,7r) with the largest suppression near (7r/2,7r/2). 

As shown in an earlier paper—, it is likely that for 
large values of U the correlations responsible for the sup- 
pressed weight in the normal state are of short-range 
antiferromagnetic origin, although the precise nature of 
the magnetic fluctuations seems less important than the 
fact that the system is close to a Mott insulator^° and 
that there are short-range spatial magnetic correlations 
induced by J. This is the stron g-co upling mechanism 
for pseudogap discussed in Refs. |43 and IH. By con- 
trast, the weak-coupling mechanism does not need to be 
close to a Mott insulator. It involves long wave length 
fluctuation3^>2i>26 that can be explained by the Two- 
Particle Self-Consistent approach. The latter mechanism 
seems to explain^^ experimental results in slightly under- 
doped Nd2-a;Cea;Cu04. For a discussion of the pseudo- 
gap at large and small U see Refs. Ilsllssl . 

In the EDC for the hole-doped case on the lower part 
of Fig. [HI one can see in the (0, 0) to (tt, tt) direc- 
tion, a gap-like feature between —t and — 2t that looks 



like the so-called waterfall feature observed in ARPES 
experiment a^^i^^'^^'^°i^^'^^ . This has some similarity to 
what has been obtained in the t ~ J model^'^ but one 
should be careful since at these energy scales it is likely 
that the one-band Hubbard model differs from more re- 
alistic multi-band models^. The gap-like feature here oc- 
curs between the spin-fluctuation induced band near the 
Fermi energy^ and the rest of the lower-Hubbard band. 
There are differences with DMFT calculations since these 
are appropriate only for higher dimension^^. 

In the hole doped case of Fig. [51 we do not have the 
resolution to address the question of whether there are 
two gaps, one associated with superconductivity, and the 
other one associated with the pseudogap, as discussed re- 
cently in the literatur o^^i^^'^^i^^ . This question is clearly 
most difficult to address in the hole-doped compounds 
where the pseudogap and superconducting gap occur in 
the same region of the Brillouin zone. In the electron- 
doped case of Fig. [3 the position of the pseudogap or of 
the antiferromagnetic gap near (7r/2, 7r/2) differs from the 
location of the maximum superconducting gap. Hence, 
an earlier VGA study'^" with better momentum resolu- 
tion allowed one to see the superconducting gap in the 
electron-doped case near (tt, 0) at the same time as the 
(pseudo)gap near (7r/2,7r/2). In our case also, the EDC 
in the lower part of Fig. [7] does have both the normal 
state gap near (7r/2,7r/2) and the d-wave superconduct- 
ing gap near (tt, 0). The MDC color scheme in the upper 
part of the same figure is on a relative scale. Hence, 
despite the fact that there is both a pseudogap and a 
superconducting gap, the MDC emphasizes that there is 
less weight in the pseudogap region near (7r/2,7r/2) than 
near the superconducting gap region near (tt, 0) , (0, tt) . 
Nevertheless, close inspection reveals that regions of con- 
stant intensity carry less weight near (tt, 0) than the cor- 
responding normal state results in Ref. [2l|. Also, the 
weight increases as one moves slightly away from the zone 
edge along the Fermi surface, refiecting the supercon- 
ducting gap. Experimentally, the superconducting gap 
in the electron-doped compounds has been seen indepen- 
dently of the pseudogap2S. 



IV. CONCLUSION 

A. Summary 

We have shown that many of the non-BCS features ex- 
hibited by superconductivity in the underdoped cuprates 
can be reproduced by CDMFT calculations for the two- 
dimensional Hubbard model. For example, the fact that 
the zero-temperature d-wave superconducting order pa- 
rameter decreases as we move towards half-filling while 
the gap in the single-particle density of states increases. 
We have demonstrated semi-quantitative agreement with 
other cluster methods and with experimental results for 
other quantities as well, such as the ARPES spectrum 
and the doping range where superconductivity is stable. 
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We have also calculated the phase diagram that de- 
scribes the competition between antiferromagnetism and 
superconductivitity for t' ~ — 0.3i for both hole and elec- 
tron doping at various values of U . There is in gen- 
eral homogeneous coexistence (superconducting antifer- 
romagnetism or equivalently antiferromagnetic supercon- 
ductivity) between antiferromagnetism and d-wave su- 
perconductivity in the underdoped region. This feature 
is quite generally observed in quantum-cluster methods, 
variational approaches and mean-field theories that do 
not allow for spin or charge density modulations on large 
length scales. In experiment, such coexistence seems to 
be the exception rather than the norm. At large U, the 
transition between the coexistence phase and the pure 
d-wave state is first order. 

Without attempting any precise fits to experimen- 
tal data, it seems that the phase diagram obtained for 
t' = — 0.3t at intermediate coupling U = 8t corresponds 
with the experimental phase diagram better than either 
smaller (U = 6t) or larger ([/ = 12t) values of U. Inter- 
estingly, adding the third neighbor hopping t" — 0.2t, 
as suggested by band structure, leads to an improve- 
ment of the qualitative agreement between the calcu- 
lated and observed phase diagrams, namely a reduced 
tendency to coexistence and a larger electron-hole asym- 
metry. The sensitivity to band parameters, expected 
at intermediate coupling where both kinetic and poten- 
tial contributions are comparable, hints towards quanti- 
tative deviations from the universal Tc {5) /T™"^ curve 
for compounds with different band parameters, such as 
La2_:rSr^Cu04 (LSCO) and YBaaCuaOy-x (YBCO). 



B. Discussion 

The major questions left open by the present work are: 
a) whether the Hubbard model by itself, without addi- 
tional interactions, additional bands, or without extrinsic 
disorder, can lead to a closer agreement with the exper- 
imental phase diagram when additional broken symme- 
tries are allowed in the space of solution s^°"'"'^°^'^°^'''^°^ b) 
if the additional broken symmetries that are present in 
cuprates appear as secondary instabilities or are essential 
to the physics of high temperature superconductivity. 

Nevertheless, the Hubbard model seems sufficient to 
lead to the overall phase diagram of the high-temperature 
superconductors and to explain many of the observed 
non-BCS features, even in the absence of a final answer 
to the last questions raised. In BCS theory, it is use- 
ful to think about the mechanism for superconductivity, 
i.e. what perturbation makes the normal state unsta- 
ble towards a superconducting ground state. Similarly, 
in the weak-coupling Hubbard model, finite tempera- 
ture BCA^^22^ and Two-Particle Self-Consistent cal- 
culations (TPSC)^, suggest that in this case it makes 
sense to think of d-wave superconductivity as being me- 
diated by antiferromagnetic fluctuation a^°^'^°^'^°^'^°^ , a 
generalization of the Kohn-Luttinger mechanism. In this 



weak- to intermediate-coupling case, pairing occurs when 
the antiferromagnetic correlation length is large22il2S. At 
strong coupling though, we have shown, in agreement 
with other approaches, that the order parameter scales 
with the superexchange interaction J, even when we do 
not allow for large antiferromagnetic correlation lengths 
(which in the cluster treatements are mimmicked by an- 
tiferromagnetic long range order). This shows that at 
strong coupling, d-wave superconductivity can occur di- 
rectly from the superexchang o^^i'^^ J, without the need 
for an antiferromagnetic "glue" or antiferromagnon ex- 
change, as advocated by AndersoniS^. Still, at large U, 
the frequency dependence of the order parameter has a 
complicated frequency dependence with multiple scales 
presentJiS. 

Antiferromagnetism and d-wave superconductivity are 
two possible phases of the strong-coupling Hubbard 
model. The normal state found in CDMFT is unstable 
to either phase or to a coexistence phase, depending on 
parameters and on doping. The CDMFT normal state 
contains both types of fluctuations, d-wave and antifer- 
romagnetic. At strong coupling it does not appear neces- 
sary to think that antiferromagnetic fluctuations lead to 
d-wave superconductivity or that d-wave superconduct- 
ing fluctuations lead to antiferromagnetism. This is per- 
haps best illustrated by the layered organic conductors 
that can be modeled by a one-band Hubbard model on 
the anisotropic triangular lattice at half-filling. There, 
one finds, theoretically and experimentally, a first or- 
der transition between antiferromagnetism and d-wave 
superconductivity^'^, so both phases are different insta- 
bilities of the same Hamiltonian. 

CDMFT and a number of methods show that maxi- 
mum pairing (maximum value of the d-wave order param- 
eter) occurs at intermediate coupling, a range appropri- 
ate for the high-temperature superconductors. So both 
weak and strong-coupling features may appeariS2iZ2ii21. 

One of the ways to distinguish weak and strong cou- 
pling limits is by the behavior of the d-wave order param- 
eter as one approaches half-filling. At strong coupling, 
it decreases as n — > 1 when U is larger than the criti- 
cal U for Mott localization. The phase fluctuations and 
short-range antiferromagnetic correlations that lead to 
the decrease need not be long range: They occur within 
the cluster. At weak coupling on the other hand, the 
decrease in the order parameter near half-filling does not 
occur unless antiferromagnetic fluctuations with correla- 
tion lengths large enough to create the pseudogap are 
included^iiiS^. Exactly at half-filling, weak and strong 
coupling limit are separated in the normal state by the 
Mott transitioniSiiii, in the absence of magnetic long 
range order. Away from half-filling, one can postulate 
that when the correlation length of the antiferromagnetic 
fluctuations become of the order of the lattice spacing, 
then J by itself suffices to lead to pairing at strong cou- 
pling. There may not be a sharp phase transition between 
the two extreme limits. 

Quantum cluster methods are powerful numerical ways 
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of attacking the problem of high-temperature supercon- 
ductivity. They clearly show that many aspects of the 
physics of the high temperature superconductors are con- 
tained in the Hubbard model. The convergence between 
the results of these methods and those of other recent 
numerical approaches suggests that we are closing in on 
methods that can provide a quantitative solution of this 
problem. 

Acknowledgements 

We thank P. Fournier and D.J. Scalapino for useful 
discussions. The present work was supported by NSERC 



(Canada), FQRNT (Quebec), CFI (Canada), CIFAR, 
the Tier I Canada Research chair Program (A.-M.S.T.) 
the NSF under grant DMR 0528969 (G.K.) and MIUR 
PRIN Prot. 200522492 (M.C.). Computations were car- 
ried out on the RUPC cluster at Rutgers, on the Elix2 
200-cpu Beowulf cluster and on the 868-cpu Dell clus- 
ter of the Reseau quebecois de calcul de haute perfor- 
mance (RQCHP), both at Sherbrooke. G.K. and A.- 
M.S.T. thank the Aspen Center for Physics where the 
final version of this paper was drafted. 



1 M. Sutherland, D. G. Hawthorn, R. W. Hill, F. Ron- 
ning, S. Wakimoto, H. Zhang, C. Proust, E. Boaknin, 
C. Lupien, L. Taillefer, et al., Phys. Rev. B 67, 174520 
(2003). 

^ M. Meinders, H. Eskes, and G. Sawatzky, Phys. Rev. B 

48, 3916 (1993). 
^ A. Macridin, T. A. Maier, M. S. Jarrell, and G. Sawatzky, 

Phys. Rev. B 71, 134527 (2005). 
* B. J. Powell and R. H. McKenzie, arXiv:0607078 (2007). 
^ G. Kotliar, S. Savrasov, G. Palsson, and G. Biroli, Phys. 

Rev. Lett. 87, 186401 (2001). 
^ B. Kyung and A.-M. S. Tremblay, Phys. Rev. Lett. 97, 

046402 (pages 4) (2006). 
^ P. Sahebsara and D. Senechal, Phys. Rev. Lett. 97, 

257004 (2006). 

T. Maier, M. Jarrell, T. Pruschke, and M. H. Hettler, 
Reviews of Modern Physics 77, 1027 (2005). 
^ A. M. S. Tremblay, B. Kyung, and D. Senechal, Low 

Temp. Phys. 32, 424 (2006). 
1° G. Kothar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. ParcoUet, and C. A. Marianetti, Reviews of Modern 
Physics 78, 865 (pages 87) (2006). 

A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

12 M. Jarrell, Phys. Rev. Lett. 69, 168 (1992). 

1^ M. H. Hettler, A. N. Tahvildar-Zadeh, M. Jarrell, T. Pr- 
uschke, and H. Krishnamurthy, Phys. Rev. B 58, R7475 
(1998). 

" M. Potthoff, M. Aichhorn, and C. Dahnken, Phys. Rev. 
Lett. 91, 206402 (2003). 

M. Potthoff, Eur. Phys. J. B (France) 32, 429 (2003). 
1® G. Biroli and G. Kothar, Phys. Rev. B 71, 037102 (2005). 
" C. J. Bolech, S. S. Kancharla, and G. Kothar, Phys. Rev. 

B 67, 075110 (2003). 
1* M. Capone, M. CiveUi, S. S. Kancharla, C. CasteUani, 

and G. Kothar, Phys. Rev. B 69, 195105 (2004). 

O. ParcoUet, G. Biroh, and G. Kothar, Phys. Rev. Lett. 

92, 226402 (pages 4) (2004). 
2° M. Civelh, M. Capone, S. S. Kancharla, O. ParcoUet, and 

G. Kotliar, Phys. Rev. Lett. 95, 106402 (pages 4) (2005). 

21 B. Kyung, S. S. Kancharla, D. Senechal, A.-M. S. Trem- 
blay M. CiveUi, and G. Kothar, Phys. Rev. B 73, 165114 
(pages 6) (2006). 

22 A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 



Phys. 75, 473 (2003). 
2^ M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 
(1994). 

2'! D. Poilblanc and D. Scalapino, Phys. Rev. B 66, 052513 
(2002). 

2^ M. CiveUi, PhD thesis, Rutgers University, (unpubhshed) 

arXiv:0710.2802 (2006). 
2^ B. Kyung, G. Kothar, and A.-M. S. Tremblay, Phys. Rev. 

B 73, 205106 (pages 13) (2006). 
2^" B. Kyung, J.-S. Landry, and A. M. S. Tremblay, Phys. 

Rev. B 68, 174502 (2003). 
2** J. Reiss, D. Rohe, and W. Metzner, Phys. Rev. B 75, 

075110 (2007). 

23 T. Maier, M. JarreU, T. Pruschke, and J. Keller, Phys. 

Rev. Lett. 85, 1524 (2000). 
^° D. Senechal, P.-L. Lavertu, M.-A. Marois, and A.-M. S. 

Tremblay, Phys. Rev. Lett. 94, 156404 (2005). 
^1 A. Paramekanti, M. Randeria, and N. Trivedi, Phys. Rev. 

B 70, 54504 (2004). 
^2 S. Sorella, G. B. Martins, F. Becca, C. Gazza, L. Capri- 

otti, A. Parola, and E. Dagotto, Phys. Rev. Lett. 88, 

117002 (2002). 

B. Edegger, V. Muthukumar, and C. Gros, Advances in 

Physics 56, 927 (2007). ^ 

P. Lee, Report of Progress in Physics l|arXiv:0708.2115P 
(2007). 

G. Kothar and J. Liu, Phys. Rev. Lett. 61, 1784 (1988). 
^® M. Inui, S. Doniach, P. J. Hirschfeld, and A. E. Rucken- 
stein, Phys. Rev. B 37, R2320 (1988). 
M. Capone, M. Fabrizio, C. CasteUani, and E. Tosatti, 
Science 296, 2364 (2002). 

M. Capone, M. Fabrizio, C. CasteUani, and E. Tosatti, 
Phys. Rev. Lett. 93, 047001 (pages 4) (2004). 
^3 R. Ofer, G. Bazalitsky, A. Kanigel, A. Keren, A. Auer- 
bach, J. S. Lord, and A. Amato, Phys. Rev. B 74, 220508 
(pages 4) (2006). 

L. Spanu, M. Lugas, F. Becca, and S. Sorella, 
arXiv:0709.2850 (2007). 
"1 S. R. White and D. J. Scalapino, Phys. Rev. B 60, R753 
(1999). 

*2 E. Pavarini, L Dasgupta, T. Saha-Dasgupta, O. Jepsen, 
and O. K. Andersen, Phys. Rev. Lett. 87, 047003 (2001). 
D. Senechal and A.-M. S. Tremblay, Phys. Rev. Lett. 92, 
126401 (2004). 



12 



B. Kyung, V. Hankevych, A.-M. Dare, and A.-M. S. Trem- 
blay, Phys. Rev. Lett. 93, 147004 (2004). 
R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, C. D. 
Frost, T. E. Mason, S.-W. Cheong, and Z. Fisk, Phys. 
Rev. Lett. 86, 5377 (2001). 

J.-Y. P. Delannoy, M. J. P. Gingras, P. C. W. Holdsworth, 

and A.-M. S. Tremblay, (unpublished) (2007). 

See Fig. 49 of Ref. i (2007) . 

B. Kyung, Phys. Rev. B 62, 9083 (2000). 

H. Mukuda, M. Abe, Y. Araki, Y. Kitaoka, K. Tokiwa, 

T. Watanabe, A. lyo, H. Kito, and Y. Tanaka, Phys. Rev. 

Lett. 96, 087001 (pages 4) (2006). 

P. Dai, H. J. Kang, H. A. Mook, M. Matsuura, J. W. 
Lynn, Y. Kurita, S. Komiya, and Y. Ando, Phys. Rev. B 
71, 100502 (pages 4) (2005). 

R. J. Birgeneau, C. Stock, T. J. M., and K. Yamada, J. 
Phys. Soc. Japan 75, 111003 (2006). 

B. Khaykovich, Y. S. Lee, R. W. Erwin, S.-H. Lee, 
S. Wakimoto, K. J. Thomas, M. A. Kastner, and R. J. 
Birgeneau, Phys. Rev. B 66, 014528 (2002). 
Y. Zhang, E. Dernier, and S. Sachdev, Phys. Rev. B 66, 
094501 (2002). 

B. Lake, G. Aepph, K. N. Clausen, D. F. McMorrow, 
K. Lefmann, N. E. Hussey, N. Mangkorntong, M. Nohara, 
H. Takagi, T. E. Mason, et al.. Science 291, 1759 (2001). 
J. Chang, A. P. Schnyder, R. Gilardi, H. M. R. nnow, 
S. Pailhes, N. B. Christensen, C. Niedermayer, D. F. Mc- 
Morrow, A. Hiess, A. Stunault, et al., Phys. Rev. Lett. 
98, 077004 (pages 4) (2007). 
J. Sonier, Rep. Prog. Phys. 70, 1717 (2007). 
E. Dagotto, Science 309, 257 (2005). 
G. T. and C. Lhuillier, Phys. Rev. B 47, 2775 (1993). 
A. Himeda and M. Ogata, Phys. Rev. B 60, R9935 (1999). 
M. Randeria, Private communication (2007). 
W.-Q. Chen, K.-Y. Yang, T. M. Rice, and F. C. Zhang, 
arXiv:0706.3556 (2007). 

A. I. Lichtenstein and M. I. Katsnelson, Phys. Rev. B 62, 
R9283 (2000). 

M. Capone and G. Kothar, Phys. Rev. B 74, 054513 
(2006). 

M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, 
Phys. Rev. B 74, 024508 (2006). 

M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, 
Phys. Rev. B 74, 235117 (2006). 

M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, 
arXiv:0707.3557vl (2007). 

M. Aichhorn and E. Arrigoni, Europhys. Lett. 72, 117 
(2005). 

T. A. Maier, M. Jarrell, T. C. Schuhhess, P. R. C. Kent, 
and J. B. White, Phys. Rev. Lett. 95, 237001 (pages 4) 

(2005) . 

T. A. Maier, M. Jarrell, and D. J. Scalapino, Phys. Rev. 
B 75, 134519 (pages 5) (2007). 

T. A. Maier, M. Jarrell, and D. J. Scalapino, Phys. Rev. 
B 74, 094513 (pages 6) (2006). 

T. A. Maier, M. S. Jarrell, and D. J. Scalapino, Phys. 
Rev. Lett. 96, 047005 (pages 4) (2006). 
A. Macridin, M. Jarrell, T. Maier, P. R. C. Kent, and 
E. D'Azevedo, Phys. Rev. Lett. 97, 036401 (pages 4) 

(2006) . 

S. Roy (2007). 

Y. Vilk and A.-M. Tremblay, J. Phys I (France) 7, 1309 
(1997). 

R. S. Markiewicz, S. Sahrakorpi, M. Lindroos, H. Lin, and 



A. Bansil, Phys. Rev. B 72, 054519 (pages 13) (2005). 
P. W. Leung, B. O. Wells, and R. J. Gooding, Phys. Rev. 
B 56, 6320 (1997). 

M. Guillot, MSc, Universite de Sherbrooke (unpublished) 
(2007). 

M. Presland, J. Tallon, R. Buckley, R. Liu, and N. Flower, 
Physica C 176, 95 (1991). 

S. D. Obertelli, J. R. Cooper, and J. L. Tallon, Phys. Rev. 
B 46, 14928 (1992). 

Y. Ando, Y. Hanaki, S. Ono, T. Murayama, K. Segawa, 
N. Miyamoto, and S. Komiya, Phys. Rev. B 61, R14956 
(2000). 

T. Kondo, T. Takeuchi, U. Mizutani, T. Yokoya, S. Tsuda, 
and S. Shin, Phys. Rev. B 72, 024533 (pages 9) (2005). 
Y. Okada, H. Ikuta, T. Kondo, and U. Mizutani, Physica 
C 426, 386 (2005). 

C. Renner, B. Revaz, J.-Y. Genoud, K. Kadowaki, and 
O. Fischer, Phys. Rev. Lett. 80, 149 (1998). 
T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). 
V. Hankevych, B. Kyung, A.-M. Dare, D. Senechal, and 

A. -M. Tremblay, J. Phys. Chem. Solids 67, 189 (2006). 

E. Motoyama, G. Yu, L Vishik, O. Vajk, P. Mang, and 
M. Greven, Nature 445, 186 (2007). 

B. P. Xie, K. Yang, D. W. Shen, J. F. Zhao, H. W. On, 
J. Wei, S. Y. Gu, M. Arita, S. Qiao, H. Namatame, et al., 
Phys. Rev. Lett. 98, 147001 (pages 4) (2007). 

J. Graf, G.-H. Gweon, K. McElroy, S. Y. Zhou, 

C. Jozwiak, E. Rotenberg, A. Bill, T. Sasagawa, H. Eisaki, 
S. Uchida, et al., Phys. Rev. Lett. 98, 067004 (pages 4) 
(2007). 

J. Graf, G. H. Gweon, and A. Lanzara, Physica C- 
Superconductivity and its applications 460, 194 (2007). 
Z.-H. Pan, P. Richard, A. V. Fedorov, T. Kondo, 
T. Takeuchi, S. L. Li, P. Dai, G. D. Gu, W. Ku, Z. Wang, 
et al., arXiv:cond-mat/0610442 (unpublished) (2006). 
W. Meevasana, X. J. Zhou, S. Sahrakorpi, W. S. Lee, 
W. L. Yang, K. Tanaka, N. Mannella, T. Yoshida, D. H. 
Lu, Y. L. Chen, et al., Phys. Rev. B 75, 174506 (pages 7) 
(2007). 

T. Valla, T. E. Kidd, W.-G. Yin, G. D. Gu, P. D. Johnson, 
Z.-H. Pan, and A. V. Fedorov, Phys. Rev. Lett. 98, 167003 
(pages 4) (2007). 

F. Tan, Y. Wan, and Q.-H. Wang, Phys. Rev. B 76, 
054505 (pages 5) (2007). 

K. Byczuk, M. Kollar, K. Held, Y. F. Yang, L A. 
Nekrasov, T. Pruschke, and D. VoUhardt, Nature Physics 
3, 168 (2007). 

M. Le Tacon, A. Sacuto, A. Georges, G. Kotliar, Y. Gal- 
lais, D. Colson, and A. Forget, Nature Physics 2, 537 

(2006) . 

K. Tanaka, W. S. Lee, D. H. Lu, A. Fujimori, T. Fujii, 
Risdiana, L Terasaki, D. J. Scalapino, T. P. Devereaux, 
Z. Hussain, et al., Science 314, 1910 (2006). 
M. Aichhorn, E. Arrigoni, Z. B. Huang, and W. Hanke, 
arXiv:cond-mat/0702391 (2007). 

M. Civelh, M. Capone, A. Georges, K. Haule, O. Par- 
collet, T. Stanescu, and G. Kotliar, cond-mat/0704.1486 

(2007) . 

N. P. Armitage, D. H. Lu, D. L. Feng, C. Kim, A. Dam- 
ascelh, K. M. Shen, F. Ronning, Z.-X. Shen, Y. Onose, 
Y. Taguchi, et al., Phys. Rev. Lett. 86, 1126 (2001). 
Q. Li, M. Hiicker, G. D. Gu, A. M. Tsvelik, and J. M. 
Tranquada, arXiv:0703357 (2007). 

E. Berg, E. Fradkin, E.-A. Kim, S. A. Kivelson, 



13 



V. Ogancsyan, J. M. Tranquada, and S. C. Zhang, 

arXiv:0704.1240 (2007). 

B. Fauque, Y. Sidis, V. Hinkov, S. Pailhcs, C. T. Lin, 
X. Chaud, and P. Bourges, Phys. Rev. Lett. 96, 197001 
(pages 4) (2006). 

S. A. Kivelson, 1. P. Bindloss, E. Fradkin, V. Oganesyan, 
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 
Mod. Phys. 75, 1201 (2003). 

M. Bcal-Monod, C. Bourbonnais, and V. Emery, Phys. 
Rev. B 34, 7716 (1986). 

K. Miyake, S. Schmitt-Rink, and C. Varma, Phys. Rev. B 
34, 6554 (1986). 



D. Scalapino, J. Loh, E., and J. Hirsch, Phys. Rev. B 34, 

8190 (1986). 

J. Carbotte, E. Schachinger, and D. Basov, Nature 401, 
354 (1999), ISSN 0028-0836. 

S. Hassan, B. Davoudi, B. Kyung, and A.-M. Tremblay, 
arXiv:0711.0739 (2007). 
P. W. Anderson, Science 316, 1705 (2007). 
"° K. Haule and G. Kothar, arXiv:0709.0019 (2007). 

A. H. Nevidomskyy, C. Scheiber, D. Senechal, and A.- 
M. S. Tremblay, arXiv:0711.0214 (2007). 



